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A new scheme for calculating masses and boost-invariant wave functions of 
heavy quarkonia is developed in a light-front Hamiltonian formulation of QCD. 
Only the simplest approximate version with one flavor of quarks and an ansatz 
for the mass gap for gluons is discussed. The resulting spectra look reasonably 
good in view of the crude approximations made in the simplest version. 
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1. Motivation for the LF Hamiltonian approach to QCD 

The method for calculating masses and wave functions of heavy quarkonia 
that is reported here stems from the program of a weak-coupling expansion 
for Hamiltonians in light-front (LF) QCD 1 . The LF form of dynamics was 
discovered by Dirac 2,3 and continues to excite imagination of physicists 4 . 
Many authors have rediscovered LF dynamics. A famous example concerns 
application to hard exclusive processes 5 . Review articles provide other ref- 
erences concerning LF formalism 6 ' 7 . One reason of the great interest is that 
the field quantization on the front hyperplane leads to 7 kinematical gener- 
ators of the Poincare group, instead of only 6 kinematical generators in the 
standard form (3 momentum and 3 angular momentum operators). Another 
reason is that the vacuum problem in the LF formulation of quantum field 
theory appears intriguingly different from the standard version. The same 
two reasons propelled also the development of the renormalization group 
procedure for effective particles (RGPEP) that is the basis of the method 
discussed here 8 . But there are two more reasons. 

The first is that the LF Fock space of free bare particles can be intro- 
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duced before one constructs the concept of a quantum field operator and 
builds Hamiltonian interaction terms for the bare quanta using such field 
operators 9 . This is useful when one attempts to mathematically define a 
theory of quarks and gluons that never appear as incoming or outgoing 
particles in scattering experiments but exist inside hadrons. Proceeding in 
this order, one can regulate the interaction terms in the LF Hamiltonian in 
a boost-invariant way. The regularization is accomplished using the relative 
transverse momenta and fractions of total " plus" momentum that the bare 
particles in interaction are carrying 8 . The transverse and "plus" momenta 
are defined with respect to the direction of the front hyperplane, the lat- 
ter conventionally defined by the condition x + = x° + x 3 = in a frame 
of reference in which the front is moving along z-axis (x 3 ), extending in 
the transverse directions of coordinates x^ . The transverse momenta of the 
particles are denoted by k 1 - and their "plus" momenta by k + = k° + k 3 . 
The regulated interaction Hamiltonian for bare particles is invariant with 
respect to boosts along the z-axis and two additional boost-like transfor- 
mations that can change transverse momenta to arbitrary values. It is also 
invariant with respect to two translations in the _L directions, translation 
in the x~ direction, and rotations around the z-axis (typically directed 
along the beam, a dominant momentum transfer, or a suitable combination 
thereof depending on a scattering experiment, but in a complete theory 
the choice should not matter). Thus, the Hamiltonian has the same struc- 
ture in a large class of frames of reference (7 dimensional). Consequently, 
one does not need to construct Hamiltonian counterterms that restore boost 
symmetry when one tries to quantitatively explain the mechanism by which 
masses, spins, and other quantum numbers of hadrons are formed. Most at- 
tractively, the basic Hamiltonian has the same structure in the rest frame of 
a hadron, where the constituent picture works 10 , and in the infinite momen- 
tum frame, where the parton model works 11 . The LF Hamiltonian approach 
raises hopes for conceptual and quantitative explanation of the constituent 
and parton models in a single and complete formulation of QCD. 

The second reason is that one can take advantage of the concept of 
potentials acting at a distance between relativistic quarks and gluons, a 
feast not conceivable in the standard approach that is defined using ob- 
jects distributed on a space-like hyperplane in space-time. Every interac- 
tion between two objects located at different points of such a hyperplane 
corresponds to a dynamical effect spreading faster than light and has to 
eventually cancel out in observables (this happens in perturbative QED 
but it is not clear how it may happen in non-perturbative QCD). There- 
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fore, it is common in the standard approach to consider only action and use 
local Lagrangian densities for field variables in a path-integral formula for 
transition amplitudes. Geometrical ideas such as strings and other nonlocal 
objects in multidimensional spaces are then used to regulate and explain 
the interaction terms in the Lagrangians. An additional argument for the 
Lagrangian approach is that it can incorporate variation of the metric in 
space-time and, hopefully, illuminate the problem of connection between 
particle dynamics and gravity 13 . But if one leaves gravity aside as too weak 
to be of an immediate concern at the scale of hadronic binding mechanism, 
it is useful to observe that the LF Hamiltonian at x + — can contain 
potential terms that act between particles separated by arbitrarily large 
distances and such interactions can obey the rule that dynamical effects 
do not spread faster than light. Namely, when the bare point-like particles 
have the same transverse positions, the four-dimensional space-time inter- 
val between them is zero no matter how large is their separation in the 
direction of x~ . In fact, the LF counterpart of the Coulomb interaction 
between two particles 1 and 2 on the LF is proportional to \x± — when 
x i = x ki an d otherwise vanishes. Precisely this type of interaction leads 
to a model of confinement in a 1+1 dimensional theory 14 . It is clear that 
the LF Hamiltonians are very singular when transverse distances between 
charged point-like particles tend to zero, and the singular terms can involve 
entire functions of the x~ distances between the particles. 

Both reasons described above indicate that one needs a powerful ultra- 
violet renormalization technique for Hamiltonians in order to develop LF 
QCD (note that the Wilsonian concept of universality could help in identi- 
fying effective Hamiltonians irrespectively of many details in setting up the 
initial bare theory). A new technique has been invented 15 ' 16 and adopted in 
a general scheme of weak coupling expansion in LF QCD 1 . More recently, 
the Hamiltonian approach has been redesigned in the form of RGPEP 8 . 
The results reported below are obtained using RGPEP and an ansatz for a 
mass-gap for gluons. 

2. Binding above threshold in heavy quarkonia 

Since there is not enough room here to thoroughly explain the RGPEP 
in application to heavy quarkonia 17 in comparison to other approaches, 
only the main steps are indicated. The central puzzle is how a systematic 
treatment of QCD can produce binding of quarks above threshold. QED 
describes quantum binding only below the mass threshold. So, how can 
binding above the threshold emerge in a rclativistic quantum theory? 
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One begins from the Lagrangian for QCD with one flavor of quarks 

1 
4 

A canonical LF procedure in gauge A + = produces a Hamiltonian with 
many terms (constraint equations are solved explicitly) 5 



C = - m)i, - -F^F^ . (1) 



H can — H,p2 + H A 2 + H^pAtp + H^^y 

+ H A 3 + H A 4 + H^AAi> + H[ dAA ]2 + £f[aAA](jK0 ■ (2) 

Each of these terms is an integral of the corresponding Hamiltonian density 
over the LF hyperplane with x + = 0, H = J dx~ d 2 x ± Tt. For example, 

n v = ^7 + ~^+ m % , n A 2 = - V(^) 2 ^ , (3) 

HipAip = g$M, = jjg 2 tpj+t a tp ^Q + j 2 ^l +t<1 ip , etc. (4) 

The fields at x + = are expanded into creation and annihilation operators 
for bare quarks and gluons, the measure is [k] — dk + cPk- 1 / '(167r 3 fc + ): 



XcUkabkace lkx + XcVkad\ ac e tkx 



(5) 



crc 



^ = V l[k] t c ela kac e- ik *+t c eZal c e ik * , (6) 



c stands for color, a for spin. The bilinear terms in TL provide kinetic ener- 
gies for the bare particles and the terms with more than two fields provide 
interactions that are regulated 17 as indicated in Section 1. The regular- 
ization implies appearance of counterterms, Hct, that restore the dynam- 
ics that was cut off by the regularization. The full regulated Hamiltonian, 
H = [H can + H CT } reg provides the initial condition for RGPEP (RGPEP 
is also used to determine Hct) 8 - The main step is to replace the canonical 
operators b, d, and a, or their hcrmitcan conjugates in Eqs. (5) and (6), 
commonly denoted by q can , by unitarily equivalent operators that create 
or annihilate effective particles corresponding to the renormalization group 
parameter A, q A = U\q can so that q x = q can and dH x /dX = [T\,H\], 
where T\ = dU\/d\ f/j. Given the initial condition = [H can + HcT] reg , 

one can systematically evaluate the Hamiltonian H\ = + ds\T Sl H s ] 
in perturbation theory. H\ is equal to H but it is expressed in terms of 
operators creating and annihilating effective particles of size 1/A with re- 
spect to strong interactions. Since H\ is expressed in terms of the creation 
and annihilation operators for effective quarks and gluons instead of the 
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bare canonical ones, it contains different interaction terms, including new 
effective potentials. For A on the order of hadronic masses, the effective 
particles are expected to correspond to the constituent quarks and gluons 
that are used to describe hadrons in particle tables 10 . The perturbative pro- 
cedure for evaluating H\ is safe from genuine infrared singularities because 
the RGPEP generator T\ is designed to exclude the possibility that energy 
denominators in perturbation theory are significantly smaller than A. 

A quarkonium eigenvalue problem for the QCD Hamiltonian H\, 
(P + H\ — P ±2 )\P) = M 2 \P), is solved by first eliminating the eigenval- 
ues P + and P 1 - of three kinematical momentum operators P^ and P^ 
(these operators are also provided by RGPEP 18 ) and obtaining an eigen- 
value equation for the quarkonium mass M (the center-of-mass motion is 
eliminated from the eigenvalue problem exactly). Still, the eigenstate \P) 
is built from the virtual effective particles in the LF Fock space and carries 
four-momentum P with P~ = (M 2 + P ±2 )/P + . In terms of the effective 
quark- antiquark, quark-antiquark-gluon, and other components: 

\P) = \Q\Q\) + \Q\Q\9\) + . . . . (7) 

This expansion may converge, in distinction from the expansion of the same 
state into canonical bare-particle sectors, because interactions in H\ are 
limited to momentum transfers smaller than A by the form factors fx that 
appear in all interaction vertices in H\. The form factors are introduced 
through the generator T\ of RGPEP. In the effective-particle basis, the 
Hamiltonian H\ takes a matrix form 



T 3 + V 3 Y 

rt t 2 + v 2 



H2+3 — 



^3 + ^ansatz ^ 

F f To + Vn 



,(8) 



in which dots denote couplings with sectors with more than 3 effective 
constituents, T refers to kinetic energy terms, V to potentials, Y to emission 
of effective gluons by quarks, and 2 and 3 to the Fock components with 2 
and 3 effective particles. The arrow indicates a truncation of the system to 
sectors \Q\Q\) and \Q\Q\g\) only, which is done at the price of introducing 
an ansatz for the gluon mass gap, 



Ansatz = \ l - ^2 J ^ ' ^ 

The ansatz is so designed that when the coupling constant a (this is the 
effective coupling at some small scale A) is extrapolated to a realistically 
strong value a s , the ansatz will be removed and the true QCD interactions 
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can be recovered order-by-order in the weak coupling expansion in a. The 
gap function [i 2 is inserted in order to model the effect of all the non-abelian 
Coulomb potentials, V3, that act in sector \Q\Q\g\) and the interactions 
that produce couplings to additional sectors with more constituents (the 
dots). It is very unlikely that the first approximation in QCD should be fi 2 = 
0. But if /i 2 ^ 0, the resulting dynamics in the \QxQx) sector is described in 
the leading order in a by the eigenvalue equation Hqq x \P) = M 2 \P), where 
the effective quark-antiquark Hamiltonian has the form (qualitatively) 

H Q Qx=T 2 + V 2 + Y^-^Y. (10) 

The main point is that the gluon emission and absorption produces di- 
verging (for small q z ) terms of the form f\ ^% — j^— 2 fx, in which the 
momentum transfer q approaches zero. This happens also in the quark 
self-interaction terms. The net effect is positive, lifting the quark energy 
above threshold. In addition, the factor dependent on [i 2 becomes 1 for 
small q irrespectively of the details of the ansatz for ^ 2 . The final result 17 
is a harmonic oscillator potential that appears as a leading correction to 
the color Coulomb interaction at typical distances between the quarks (the 
Coulomb term appears with the Breit-Fcrmi spin factors). Technically, it 
is the harmonic oscillator term that leads to the binding above threshold, 
M > 2m, where m is the mass ascribed to the quarks. Such effect is ab- 
sent in positronium in QED because there are no Coulomb-like interactions 
between photons and electrons and no mass-gap for photons . 

3. Masses and wave functions in the crudest approximation 

The resulting eigenvalue equation for quarkonium wave function can be 
solved numerically and the mass spectrum depends on the choice of the 
coupling constant a and quark mass m at some value of A. The Breit-Fermi 
terms include three-dimensional <5-functions that are smeared and made 
finite by the presence of the form factors fx- If one assumes olm z ~ 0.12, 
the RGPEP evolution with one flavor of quarks in the same Hamiltonian 
scheme 19 gives a ~ 0.326 at A ~ 3.7 GeV. Table 1 shows masses of bb 
quarkonia obtained for a = 0.326 and m = 4857 MeV at A = 3699 MeV, 
adjusted to fit masses of xi(lP) and xi(2P). The pattern of differences 
in the 4th column agrees with expectations in the new scheme. All details 
concerning this calculation and results for some other quarkonia can be 
found elsewhere 20 . The oscillator frequency corresponding to Table 1 is w 
= 182 MeV. The bottom line is that the realistic value of a is near 1/3 and 
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Table 1. Qualitative illustration of bottomonium masses. 



meson name 


calculation (MeV) 


experiment (MeV) 


difference (MeV) 


T10865 


10725 


10865 


-140 


T10580 


10464 


10580 


-116 


T3S 


10382 


10355 


27 


X22P 


10276 


10269 


7 


Xi2P 


10256 


10256 





X02P 


10226 


10232 


-6 


T2S 


10012 


10023 


-11 


X2lP 


9912 


9912 


-1 


XI IP 


9893 


9893 





XolP 


9865 


9859 


5 


T1S 


9551 


9460 


91 


% IS 


9510 


9300 


210 



the new Hamiltonian approach to QCD can be further studied in a weak 
coupling expansion in the case of heavy quarkonia, including many effects 
in the complex relativistic color dynamics of virtual quarks and gluons. 
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